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Abstract 

We investigate various aspects of the integrability of the vertex models associated 
to the D% affine Lie algebra with open boundaries. We first study the solutions of the 
corresponding reflection equation compatible with the minimal symmetry of this system. 
We find three classes of general solutions, one diagonal solution and two non-diagonal 
families with a free parameter. Next we perform the Bethe ansatz analysis for some of 
the associated open D\ spin chains and we identify the boundary having quantum group 
invariance. We also discuss a new D\ i?-matrix. 
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1 Introduction 



Much work has been done in integrable lattice statistical mechanics models with open boundary 
conditions, since Sklyanin U generalized the quantum inverse scattering method to tackle the 
boundary problem. The bulk Boltzmann weights of an exactly solvable lattice system are 
usually the non-null matrix elements of a .R-matrix R(X) which satisfies the Yang-Baxter 
equation. The integrability at boundary, for a given bulk theory, is governed by the reflection 
equation, which reads 

Rn(X - aO K- (X)R 2 i(X + n) K- 0) =K- {^)Rn{X + //) K- (X)R 2 i(X - //) (1) 

where the matrix 2f_(A) describes the reflection at one of the ends of an open chain. Similar 
equation should also hold for the reflection K + (X) at the opposite boundary. However, for 
several relevant lattice models K + (X) can be directly obtained from lf_(A). For example, this 
is the case of models whose R(X) matrix satisfies extra properties such as unitarity, P and T 
invariances and crossing symmetry fl], ^ . 

Therefore, the first step toward constructing integrable models with open boundaries is 
to search for solutions of the reflection equation. To date, solutions of this equation have 
been found for a number of lattice models ranging from vertex systems based on Lie algebras 
!§, [|, ||, U to solid-on-solid models and their restriction . Classification of such solutions for 
particular systems || as well as extensions to include supersymmetric models |J can also be 
found in the literature. 

In spite of all these works, there is an interesting vertex model based on the non-exceptional 
D 2 Lie algebra for which little is known about the solution of the corresponding reflection 
equation. This is probably related to the fact that the D 2 i?-matrix does not commute for 
different values of the rapidity consequently the trivial diagonal solution i£_(A) = / does 
not hold for this system ||. The purpose of this paper is to bridge this gap, by presenting 
what we hope to be the minimal solution of the reflection equation for D 2 vertex models. This 
result offers us the possibility to understand a relevant open problem which is the integrability 
of the D 2 n vertex model with quantum algebra symmetry. In fact, this symmetry has been 
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found for all vertex models based on non-exceptional Lie algebras except for the D\ model. 
It turns out that, by earring out a Bethe ansatz analysis, we are able to identify this symmetry 
for the simplest D\ model and conjecture it for arbitrary values of n. 

We have organized this paper as follows. We start next section by considering the reflection 
equation for the D\ vertex model. We find one diagonal solution without free parameters and 
two non-diagonal families which depend on a free parameter. We also derive the corresponding 
integrable one- dimensional open spin chains. In section 3 we present the Bethe ansatz solutions 
of the open D\ spin chain associated to the diagonal i^-matrix and to a special manifold of 
the first non-diagonal family. This allows us to identify the quantum group symmetry for the 
D\ model. In section 4 we generalize the .fT-matrices results of section 2 for arbitrary values 
of n > 2. Section 5 is reserved for our conclusions as well as a discussion on possible new 
i?-matrices. In Appendix A we collect some useful relations and Appendix B contains a new 
D\ i?-matrix as well as its boundary behaviour. 



2 The D\ i^-matrices 

The D\ vertex model has four independent degrees of freedom per bond and its Boltzmann 
weights preserve only one U(l) symmetry out of two possible ones. Here we are interested in 
looking at solutions of the reflection equation that commute with this symmetry. We find that 
the most general .fT-matrix having this property is 

\ 

Y 2 (X) Y 5 (X) 
Y 6 (X) Y 3 (X) 
Y 4 (X)j 

Our next step is to substitute this ansatz in equation (|T]) and look for relations that con- 
straint the unknown elements Yj(X),j = 1, . . . , 6. Although we have many functional equations, 
a few of them are actually independent, and the most suitable ones have been collected in Ap- 
pendix A. The basic idea is to try to solve such equations algebraically, which hopefully will 



fYiW 



V o 
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produce a general ansatz for functions Yj(X) containing several arbitrary parameters. The 
general strategy we use is to separate these equations in terms of ratio of functions depending 
either on A or on fi. From the relations (A.5-A.7) one easly concludes that the simplest possible 
solution is to take Y 5 (X) = Y 6 (X) = 0. This is the diagonal solution, and by employing the 
"separation variable method" described above for the relations (A.5-A.7) we are able to fix the 
following ratios 

Y 2 (A) _ e A — Pi Y 3 (A) _ e x -p 2 Y 4 (A) _ e A - fa 

Yi(\) e~ x — Pi ' Yi(\) e~>-/V Y 3 (X) e" A - /3 3 U 

where Pj, j = 1,2,3 are arbitrary constants. 

These relations enable us to write an ansatz for three unknown functions in terms of a 
normalizing factor, say Yi(X). Substituting the relations (3) back to the reflection equation 
(|I|), we conclude that all the parameters ftj are fixed by 

Pi = -P2 = -l/Pz = ^ ( 4 ) 



where q is the deformation parameter of the D\ i?-matrix [10]. This leads us to our first 
solutions with no free parameter, 

y«(A) = 1, K<"(A) = tz\, K 3 « (A ) = ^±4, r«"(A) = ; A A + + ^ *rjf (5) 

Next we turn our search for non-diagonal solutions now with both Y${X) and Yq(X) non 
null. From the equations (A.8-A.11), we notice that it is possible to solve Y^A), 13(A) and 
le(A) in terms of 15(A). At this point we should keep in mind that we are looking for regular 
i^-matrices, i.e. K_(0) ~ identity. After some simplifications, we find the following general 
solutions 

ei(l + e 2A )[y- 5 (A) + F 6 (A)] = e 2 e A [F 5 (A) - Y 6 (X)} (6) 
(e 2A - 1)[Y 2 (A) + y 3 (A)] = e 3 e A [y 6 (A) - y 5 (A)] (7) 

y 2 (A)-y 3 (A) = e 4 [y 5 (A)-y 6 (A)] (8) 
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where 6j are four arbitrary parameters. These are linear equations which can be easily solved 
for the ratios Y 2 (\)/Y 5 (\), F 3 (A)/F 5 (A) and Y e (\)/Y 5 (\). Taking this into account as well as 
equations (A. 5) and (A. 7), we end up with the following ansatz for functions Yj(X) 



Yi(A) = (e 5 e 2A + e 6 e A + e 7 )/e A , K>(A) = (1 + e 2A ) e 4 (e 2A - 1) - e 3 e 



2A\ 



_2A 



Y 3 (X) = (1 + e 2A ) -e 4 (e 2A - 1) - e 3 e A , F 4 (A) = (e 8 e 2A + e 9 e A + e 10 )e 



,2A 



.2A 



„2A 



1) e 2 e A + e 1 (l + e 2A ) , Y 6 (X) = (e 2A - 1) e2 e A - ei (l + e 



(9) 



(10) 
(11) 



having altogether ten free parameters. Substituting this ansatz back to the reflection equation 
and after involving algebraic manipulations, we find that nine parameters are in fact fixed, 
leading us to two classes of non-diagonal solution with a free parameter. The first class is 
given by 



yi (2) (A,e_) = (e 2A + q)(eqe 2X - l)e~\ Yf>(\,£-) = (e 2A + q)(£q ~ e 2X )e 



,2A 



-2 2A 



A v (2). 



„2A 



.2A\„A 



n (2) (A,e_) 



(1 + e 2A ) 



2(e 2X -l)^q-e\l + q)(l-eq) 



(l + e 



2A\ 



r 5 (2) (A,e_)=r 6 (2) (A,e 

while the second family is 



-2(e 2A -l)^g-e A (l + g)(l-e<?) 
(e 2A - 1 



[l-qW-q + iy 



Y} 3 \\,Z-) = (e 2A - g)(£_e 2A - l)e" A , Y?\\Z-) = (e 2A - q ){£_ - r"), 



A v (3)/ 



,2A\„A 



n (3) (A,e-) 



e 2A - 1) 



^ 6 (3) (A,e- 



(e 2A - 1) 



2(e 2A + l) v /^<7 + (l + g)(l + £-y 



2(e 2A + l)J^q+(l + q)(l+^y 



(12) 
(13) 
(14) 
(15) 

(16) 
(17) 
(18) 
(19) 



where £_ is an arbitrary parameter. 
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Since the D\ i?-matrix is PT invariant and crossing symmetric, the K + (X) matrices at the 
opposite boundary are easily derived from the above solutions 0, [J. More precisely, we have 



K + (\,£ + )=K t _(\n[q]-\,Z + )M 



(20) 



where M is a matrix related to the crossing matrix V by M = V f V 0. From the results of 
Appendix A, we have that for the D\ model M is given by 



M — diag(g, 1,1, q' 1 ) 



(21) 



Having found the K±(\) matrices, one can construct the corresponding commuting transfer 
matrix r(A). Following Sklyanin M, we have 



t (i,m)( X ) = Tr a 



a (m) a (0 , 

K + (A)T(A) K_ (\)T-\-\) 



(22) 



where T(A) = R a i{X) ■ ■ ■ R a x(X) is the monodromy matrix of the associated closed chain with L 
sites. This means that the three families of if- (A) matrices we found will produce nine possible 
types of open boundary conditions. The corresponding Hamiltonian of the spin chains with 
open boundaries are obtained by expanding the transfer matrix t^ l ' m \\) in powers of A. When 
Tr[i4 m) (0)] is non-null, the Hamiltonian H^ l ' m ' is proportional to the first-order expansion PJ 



L-l 



H {l ' m) = E H k,k+i 



k=l 



a (0 

1 d K_ (A) 
2C d\ 



A=0 + 








(m) 




Tr 

±L a 


(0)H La 


Tr 







(23) 



where H^k+i — -Pfc,fc+i^-Rfc,fc+i(A)|A=o is the two-body bulk Hamiltonian and ( is the normal- 
ization R 12 (0) = CPi 2 Q 

For the first two solutions we indeed have Tr[lf+(0)] 7^ while for the third one Tr[iT + (0)] = 
0. In this last case one has to consider the second order expansion in the spectral parameter A 
[ nj| . We find convenient to write the expression for the Hamiltonians in terms of Pauli matrices 
a^i and a z ai with components a =f, I acting on the site % of a lattice of size L. In terms of 



1 The normalization we use for R(X) (see Appendix A) produces £ = (q — l/q) 2 for D\ model. 
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these operators and up to irrelevant additive const ant sQ, we have 



fc=l 



<*=T4 



- E A*S B) (e+)<L+* (m W,w^ + ^ (m) (e + )4^j; L + j{ m) (e + )< LffT - L l 



(24) 



where the expression of the bulk part Hk,k+i is 
(g - 1/g) 



H, 



k,k+l 



JrCr ],k a ],k+l a l,k (T l,k+l + 'T.fc a 'T J fc+l Cr 4.* Cr 4>fc+l 



,<J f h a 



+ 1 1 o JK,fc a 4,fc + a T,fc+i a 4,fc+i) 1 (°T,fc°4,ifc+i + °4,fc°T,fc+i> 



1 \2 



(g + l/g) 1 3 



4 

(g - i/g) 



E «* - <fc+i) - 2 (<? + 



«=T4 



Turning to the boundary interactions we found that the chemical potentials are given by 



^(0 = 1/2-1$, ^(0 = 1/2 + /^ 



/^t (0 — (l+</)(CV-l) 

^ ) K) = ^ 3) = i 



(l+g+2gg) 



(l+q-2^) 

'(i+<?)(£V-i) 



(25) 



2 We also note that we have normalized the Hamiltonian by the pure imaginary number. 
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while the on-site parameters 8® and J$p(£) are 



6 (i) 



5 m 



_ (9- 


1) 


- 2(1- 


H) 


_ (ff- 


i) 


2(1- 


H) 


_ (H 


-i) 


2(9- 


-i) 



(26) 



and 



j t (1) (0 = J 

</| 3) (0 



(i) 

j. 

K 2 )( 







(27) 



(l+9)(l+g)-4y/gi 
(9-l)(«-l) 



rW/t^ _ rWf^ _ (g-l)(l+g 2 g) 
J T W - J | KQ) ~ (l+ ff )(fa g -l) 

r( 3 )^ _ ( 1 +g)( 1 +0+ 4 v / gj t(3) /t\ 
(?-!)(?-!) ' J i ^ 

A natural question to be asked is which (if any) of these solutions would lead us to an 
integrable D\ model with quantum algebra symmetry. One way to investigate that is by 
applying the Bethe ansatz method to diagonalize the above open spin chains. This allows 
us to extract information about the eigenspectrum, which in the case of quantum algebra 
invariance, should be highly degenerated (see e.g. ||12|| ). In next section we will discuss this 
problem in details. 



3 Bethe ansatz analysis 

The purpose of this section is to study the spectrum of some of the open spin chains presented 
in section 2 by the coordinate Bethe ansatz formalism. One of our motivations is to identify 
the boundary that leads us to the quantum group symmetry. We begin by noticing that the 
total number of spins N s = Yh=i E a =t,| * s a conserved quantity and its eigenvalues ns 
labels the many possible disjoint sectors of the Hilbert space. Therefore, the wave function 
solving the eigenvalue problem H \ ^ ns ) = E"' m '(L) \^ ns ) can be written as follows 

l*™> = E E / (ai '"^W • • • , *QJ<* Ql • • • < s , XQns |0> (28) 

a j x Qj 

where |0) denotes the ferromagnetic state (all spins up) and 1 < xq 1 < xq 2 < • • • < XQ ns < L 
indicate the positions of the spins. 
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We will start our study by first considering the open spin chain if t 1 ' 1 ) corresponding to 
the diagonal J\~-matrix solution. As it is customary we begin our discussion of the eigenvalue 
problem in the sector of one down spin, ns — 1. In this sector, we find that for 1 < x < L 

2(l/ q - q ) E(1,1){L)f(a){x) = (^-2)A/ (a H^)+/ (Q) (^+l)+/ (a) (^-l)+^ 1 4^/ (a) (^) ! « =U 

(29) 

where we have defined A = q + 1/q. The matching condition at the left and right boundaries 
gives us the following constraints 



and 



' / (T) (0) 

v / (i) (o) 



K f {l HL + l) 



/ (T) (1) 



(30) 



where the matrices parameters are given by 

3I + v /g + g 2 (I + 3 v /g) 
= TTT, — — — ' Pn 



Via = 



4g(/ + Vg) 

3 - 2g + 3q 2 - 21 \fq + llq^fq 

4g 



/u>(l) ; 

-3/ + v ^ + g 2 (-/ + 3 v ^) 
4< Z (-/ + v /g) 

3 - 2g + 3g 2 + 2/^ - 2/g^ 



(31) 



4g 



dit — (iii — —djA — —dr,\ — 



q - 1/q 



(32) 



(33) 



(34) 



In order to go ahead it is crucial to notice that both boundary constraints (30) and (31) can 
be diagonalized by the same unitary transformation U. After performing this transformation 
the new components f a (x) = Uf a (x) satisfy 



and 



( 7 (T) (o) \ 
{ / a) (o) ) 

f m (L + l) 
/U>(L + 1) 



1 



^ 

u 2 







2 

1 



/ (T) (1) 
/ a) (l) , 

v f a) (L) , 



(35) 



(36) 
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Clearly, equation (29) for 1 < x < L remains the same but now for the transformed 
amplitudes f^{x). Now we reached a point in which one can try the usual Bethe ansatz (e.g. 
see ref. |13], namely 

= A a (k)e ikx - A a (-k)e- ikx (37) 
and by substituting this ansatz in (29) we obtain the following eigenvalue 

— *—Ep±\L) = (L - 2) A + 2 cos(fc) + ^£ (38) 
2(1/9 -q) 4g 

The fact that this ansatz should be also valid for the ends x — 1 and x = L provides us 
constraints for the amplitudes A(k) and A(—k), which reads 

A{-k) = -e lk A{k) and A(-k) = 2 -^-^-e 2l(L+1)k A{k) (39) 

(1 — -^e 1 ) 

whose compatibility gives a restriction on the momentum k, namely 

( P ik _ A\ 

e" fcL /f , \\ = 1 (40) 



The next task is to generalize these results for arbitrary numbers of down spins. For a 
general multiparticle state, we assume the Bethe ansatz wave function 

ns 

\x Ql , ■ ■ -,x Q J = $>n(P) -Q e i^ Qj ] A{kpQij ... } k PQNe ) aQi ,.., aQns (41) 

P j=l 

where P is the sum over all the permutations of the momenta, including the negations kj — > 
—kj, and the symbol sgn accounts for the sign of the permutations and negations. It turns 
out that for configurations such that \xQ i — xq.\ > 2 the open spin chain H^ 1 ' 1 ' behaves as a 
free theory and the corresponding eigenvalues are 

-E^ x,x \L) = (L — 1)A + ^[2 cos(kj) — A] + — — V- (42) 



2(1/9-9) ' ' " ' p{ L " JJ 1 4g 
The new ingredient for ns > 2 is that the nearest neighbor spin configurations enforce 
constraints on the amplitude of the wave function. This condition enhances a relation between 
the exchange of two states such as {(h, aj); (kj,aij)} and {(kj,aij); (ki, a^)} which ultimately 
is represented by the two-body scattering 

A... aj,ai— (' ■ ■ ? kj., ki, ■ ■ ■) = Sij(ki, kj)A... auCej ...(- ■ ■ , ki, kj, ■ ■ •) (43) 
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while the reflection at the left and right ends generalizes equation (29), which now reads 

= •^... (/•>---) (44) 

~-kj) = jf^f^e^ L+ ^A... at X---,k 3 ) (45) 
(1 — f e 1 ^) 

Fortunately, the bulk two-body scattering amplitude Si y j(ki, kj) has been recently identified 
in ref . |T5] for the periodic chain. This result is of enormous help here since it allows us to choose 
the suitable parametrization for the momenta kj in terms of the S"-matrix rapidities Xj, which 
is 



sinh(Aj + iy/2) 

where we have conveniently defined q = e 1 " 1 . For explicit expression of the non-null S'-matrix 
elements see ref. HI"5f . 

In this general case, the compatibility between the bulk and boundary scattering constraints 
(43-45) leads us to the Bethe ansatz equation for the momenta kj 

( ikj _ a\ 

e 2lk > L ( A lk . 2 > =AAki,---,k ns ) (47) 

where Aj(ki, • ■ ■ , k ns ) are the eigenvalues of the auxiliary inhomogeneous transfer matrix tj = 
Sj ns (kj, k ns ) ■ ■ ■ Sji(kj, ki)Sij(kl, —kj) ■ ■ ■ S ns j(k ns , —kj). The integrability of this latter inho- 
mogenous problem follows from the fact that the 2x2 identity iT-matrix is a solution of the 
reflection equation associated to the two-body scattering Sij. As was shown in ref. [113] there is 



no need of a second Bethe ansatz to solve this auxiliary eigenvalue problem. By adapting the 



results of ref. [J5j to our case and by relating the momenta kj and the rapidities Xj by equation 



(46) we find that the Bethe ansatz equations are given by 



sinh(Aj — iy/2) 
sinh(Aj + iy/2) 



-i 2L 



cosh(A i + iy/2) ™ sinh(A i /2 - X k /2 - iy/2) sinh(Aj/2 + X k /2 - iy/2) 

cosh(A j - iy/2) ^ sinh(A j /2 - A fc /2 + iy/2) sinh(A i /2 + A fc /2 + iy/2) 

j = l,~-,ns (48) 



and the eigenvalues (42) in terms of the rapidities Xj are 

ns i 

E^\L) = -8sin 3 ( 7 ) —T-, T7^ " < L ~ ^ sin ^ + 4 sin( 7 ) sin 2 ( 7 /2) (49) 

~^ cos(7 i ) — cosn(2Aj) 
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Clearly, the Bethe ansatz equations for the open spin chain H^ 1 ' 1 ^ are not just the "doubling" 



of the corresponding results of the closed chain with periodic boundary conditions Jl6], [OJ due 
to an additional boundary left hand factor. We recall here that the "doubling" property has 
been argued |l7j to be one of the main features of a quantum algebra invariant open spin chain 
at least for standard forms of comultiplication. Looking at the spectrum of H^ 1 ' 1 ^, however, 
we notice a certain pattern of degeneracies which suggests an underlying hidden symmetry. It 
could be that the diagonal boundary solution corresponds to an asymmetric form of coproduct 
since this, in principle, is allowed too p8|| . 

Next we turn our attention to the first non-diagonal solution and its corresponding open 
spin chain. In this case, at least for generic values of £±, the Bethe ansatz construction we 
just explained above needs further generalizations. This can be seen even at the level of 
one down spin state, since there is not a unique transformation that diagonalizes both left 
and right boundary matrix problems. However, there is a particular manifold, £ + = in 
which our previous Bethe ansatz formulation is still valid. Fortunately, as we shall see below, 
this special manifold will be sufficient to single out the boundary leading us to the quantum 
algebra symmetry. Since for £ + = the Bethe ansatz analysis is very similar to the one 
just described above, we restrict ourselves to present only the final results. We found that the 
Bethe ansatz equations for the Hamiltonian H^ 2 ' 2 ^ at £ + = g£_ are 

sinh(A j - i7/2)l _ -p- sinh(A j /2 - \ k /2 - i 7 /2) sinh(A j /2 + A fc /2 - z 7 /2) ^ 
sinh(A j +i 7 /2)J ~ 1} X sinh(A i /2 - A fc /2 + z 7 /2) sinh(A i /2 + \ k /2 + 1 7 /2) ' 3 ~ r " ,ns 

(50) 

while the corresponding eigenvalues are given by 



ns i 

E^\L) = -8sin 3 ( 7 )]T^- ^^_ 4 (L-l)sin(2 7 )-4sin( 7 ) 

~[ cos( 7 J — cosh(2Aj) 



E(^ 2) (e-)-^ 2) (e + )) + 2^ 

(51) 



Now the Bethe ansatz equations do have the "doubling" property at £ + = g£_ and this is 
an extra motivation to investigate the eigenspectrum of H^ 2,2 '. It turns out that at the value 

= and therefore £ + = we discover that the spectrum of the open chain if ( 2 > 2 ) is specially 
highly degenerated. In fact, after some algebraic manipulations, we check that for £± = 
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the Hamiltonian if ( 2 > 2 ) has the appropriate boundary coefficients to ensure commutation with 
U q (D\). Therefore, we finally managed to identify the quantum algebra symmetry for the D\ 
vertex model. 

Finally, it seems desirable to solve the open spins chains associated to the non-diagonal 
solutions for arbitrary values of the parameters £±. The coordinate Bethe ansatz method, 
however, leads us to cumbersome calculations even for the first excitation over the reference 
state. In such general case it seems wise to tackle this problem by using a more unifying 
technique such as the algebraic Bethe ansatz approach. Since the basics of this method has 
been recently developed for the D\ vertex models [fTS] we hope to return to this problem 
elsewhere. 



4 The D z n i^-matrices 

Here we shall consider the generalizations of the i\~-matrices solutions of section 2 for the 
general D\ model. This system has n — 1 distinct U(l) conserved charges, and the i\~-matrix 
ansatz compatible with these symmetries can be represented by the following block diagonal 
matrix 

K_(X) = diag(r 1 (A), • ■ ■ , F n _!(A), A(X), Y n+2 (X), ■ • • , Y 2n (X)) (52) 
where A (A) is a 2 x 2 matrix 

/ Y n (X) Y 2n+1 (X)\ 
A(X) = (53) 
V>W 2 (A) Y n+1 (X) ) 

where Yj (A), j = l,---,2n + 2 are functions we have determined by solving the reflection 
equation. Notice that for n = 2 we recover our starting ansatz of section 2. 

Substituting this ansatz into the reflection equation, we realize that the simplest possible 
solution is the symmetric one, namely 

Yi(A)=y a (A) = ... = y^ a (A) and Y n+2 (X) = Y n+1 (X) = ■ ■ ■ = Y 2n (X) (54) 

It turns out that the remaining functional equations for the functions ^(A), Y n (X), Y n+ i(X), 
Y 2n (X), Y 2n+1 (X) and Y 2n+2 (X) are very similar to those presented in the appendix A. Therefore, 
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they can be solved by the same procedure described in section 2 and in what follows we only 
quote our final results. As before we find three general families of fT-matrices, and the diagonal 
one is given by 



y} 1 \\) = i, yW(A) 



e A _ /g -(n-l)/2 
e -A _ J q -{n-l)/2 



e A + /g-(n-l)/2 (1) e A + /g-(n-l)/2 g A _ /g (n-l)/2 

!> *2n (A) - "I 



e -A + J 5 -(n-l)/2 ' * W e -A + J 5 -(n-l)/2 g-A _ J q (n-l)/2 

The one-parameter families of non-diagonal J\~-matrices are given by 



(55) 
(56) 



n (2) (A,£- 



r 2 M ry" l )(£q»-\-™ -l)f-\ Y^(X.L 

e 



'l 1 ^2A^ 



*2n+l(A, £-) - *2n+2(A, £-) — 



2 

;i + e 2A ) 

2 

(2) 



2 e 



e 2A + g n-l ) ^2 g n-l_ e 2A )e A ^ 
2A , , ( A,1 . f/ " •),;, ^ i ; 

2A i\t „n-l A/ 



-2(e 2A - l)^"" 1 - e A (l + ^(l - e-q"- 1 ) 

2A_!) 



1-5 +l)e 



A 



and 



n (3) (A,e-) = (e 2A - (T 1 )^ - l)e- A , (A,e-) = (e 2A - g"- 1 )^- - e 2A )e 



„2A 



,2A 



(3), 



„2A 



,2A\ A 



^2n+l(A, C-) 



>; {3) (A,e-) = v; ( ?i(A ) e-) = 

(e 2A - 1) 



(l + e 



2A\ 



l-9 B_1 )(e--l)e' 



^2n+2(A, £-) — 



2 

; 2A _ ^ 



(58) 
(59) 
(60) 

(61) 
(62) 

(63) 

(64) 



2(e 2A + l)^_g (n - 1)/2 + (1 + ^^(l + £-)e A 

-2(e 2A + 1) ^fi~-q {n - 1)/2 + (1 + g n " 1 )(l + i-)e x 

The next natural step is to search for asymmetric J\~-matrices for n > 3, i.e. those having 
Y 1 (X) ^ Y 2 (X) ± ■ ■ -y n _i(A) and F n+2 (A) ^ F n+3 (A) ^ • • •Y 2n (X). In this case the number of 
free parameters grows rapidly with n and the solution of the reflection equation becomes more 
involving. To illustrate that, we consider the _D| model and for sake of simplicity we look first 
for diagonal solutions. There are six functions Yj(X) to be determined and their ratios are fixed 
by choosing some easy looking relations coming from the reflection relation. More precisely, 
we have found the following equations 

y„m p a _ rn vJW P A _ ro 

(65) 



Y 2 (X) e 2X - Cl Y 3 (X) e x -c 2 F 4 (A) e A - c 3 



Yi(A) e- 2A - Cl ' Y 1 (X) 



C2 n(A) e- 



c 3 



13 



y 6 (A) _ e A - c 4 F 5 (A) _ e- 2A - c 5 

y 4 (A) e- A -c 4 ' y 6 (A) e 2A -c 5 1 J 

where Cj are once again constants yet to be determined. Substituting these relations back to 

the reflection equation we find only one possible manifold for the parameters Cj, which reads 

c\ = C5 = —1 and C2 = C3 = C4 = (67) 

After an appropriate normalization, this solution leads us to a new diagonal f^-matrix for 
the D\ model 

ir Di (A) = diag(e- 2A ,l,l,l,l,e 2A ) (68) 

It is plausible that this "almost unity" solution and its extensions generalizes for arbitrary 
values of n > 4. Next we have looked at the possibility of asymmetric non-diagonal solutions 
for the -D 2 model. It turns out that, within our algebraic approach, we did not found any of 
such solutions. However, this possibility should not be completely rule out, at least for general 
n, since we have so many free parameters that the chance to miss a particular integrable 
manifold is high. In general, classification of the solutions of the reflection equation seems to 
be an intricated problem even for simpler models ||. We hope, however, that our X-matrices 
results prompt further investigation concerning this problem for the D 2 vertex models. 

We would like to conclude this section with the following remarks. The K + (X) matrices 
can be obtained from K-(\) by the isomorphism 

K + (X) — Kt[(n — 1) ln[g] — X]M (69) 

where M is a diagonal matrix given by 

M = diag(g^ 3 ), q 2n ~\ ■ ■ ■ , 1, 1, • • • , q-^\ q ~^) (70) 

Once we are equipped with K±(\) matrices, the construction of the corresponding open 
spin chains is possible along the lines of section 2. Similarly, at least for the diagonal solution, 
one can also repeat our Bethe ansatz construction without further technical difficulties. In 
particular, we conjecture that the open spin chain associated to the first non-diagonal solution 
at £± = is the one having the underlying quantum group symmetry. 
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5 Concluding Remarks 



In this paper we have made a great deal of progress towards the understanding of the integra- 
bility of the D 2 n vertex model with open boundaries. We have investigated the solutions of the 
associated reflection equation and found three general families of i^-matrices which respect 
the minimal U(l) symmetries of this system. We have carried out a Bethe ansatz analysis for 
the simplest case, D\ model, revealing to us that the first non-diagonal solution at £± = 
possesses the special quantum algebra symmetry. In fact, the structure of the ^-matrices at 
this particular point leads us to conjecture that this will be the case for arbitrary values of n. 

We believe that our results open an enormous avenue for further investigations. One clear 
possibility is to use the Bethe ansatz results of section 3 to compute the thermodynamic 
behaviour, the bulk and the surface critical exponents. It would be also interesting to generalize 
our results of section 3 for all sort of open boundary conditions and for arbitrary values of n. In 
this case, probably the most suitable tool would be instead the algebraic Bethe ansatz approach. 
This method would allow us to show that indeed the Bethe ansatz states are highest weight 
states of the underlying quantum algebra in the case of the first non-diagonal family at £± = 0. 

Other interesting issue is to apply the notion of the quantum group twisting [19| to find 
out slightly different D\ R- matrices. As a result, this might lead us to integrable models with 



very different behaviour, for an example see ref . |20j . The practical implementation of twisting, 
however, seems to be quite involving specially for an algebra such D 2 n . To shed some light to this 
problem we proceed in a much more phenomenological way. Motived by the structure of the 
non-diagonal solutions, we add extra Boltzmann weights to the Jimbo's i?-matrix to account 
for such boundary terms at the level of the associated bulk Hamiltonian. Next step is to try to 
solve the Yang-Baxter equation for this novel i?-matrix structure. It turns out that we succeed 
to find a new i?-matrix solution for the D\ model. Since this involves many technicalities, we 
have summarized it in appendix B together with the study of the corresponding solutions of 
the reflection equation. We hope that these results will be useful to motivate further progress 
in this problem. 
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Appendix A : R matrix properties and reflection equations 
In this appendix we briefly discuss some useful properties of the D 2 i?-matrix [[K]]. We 
also present for n = 2 some relevant relations derived from the reflection equation. The D 2 n R- 
matrix satisfies, besides the unitarity and regularity, extra relations denominated PT invariance 
and crossing symmetry, PT-Symmetry: 

Pi2Ri2(X)Pi2 = R\ 1 2 t2 (X) (A.l) 

Crossing-symmetry: 

R ™W = TVi rfr r i xi V R ^ n - ^ - A ) V 1 (A.2) 

Cl(n - l)ln[gj - AJ 

where £(A) is a normalization function and V is the following crossing matrix 

V = antidiag |g-( 2 «- 3 )/ 2 , q ~^)/\ . . . ; i_, 1, 1, ^ . . . , g (a«-6)/2 > g (2n-3)/2 j (A 3) 

Here we find convenient to normalize the original Jimbo's i?-matrix by an overall factor 
e 2X q n and the function £(A) is given by 

A (n-l) 

C (A) = (^-e- A )(^-^ r ) (A.4) 

Next we present the simplest relations derived from the reflection equation we used in 
section 2. For sake of simplicity we shall use the following notation Y{(x) = Yi,Yi(y) = 
Y i ,Wj(x — y) = Wj,Wj(x + y) = w'j. Considering this notation, the relations we have selected 
from the reflection equation are given by 



Wn 



-wsY^ + w±Y x Y 2 - w h Y x Y h + w 6 YiY 6 
16 



Wo 



Wo 



—W 



wXYi + w' 5 (Y;Y 5 + Y 2 Y 6 ') + w' 3 (Y 2 % + Y 5 Y,' 
' -w 3 YlY 3 + w 4 Y x Y 3 + wgYX - w 5 Y;Y e : 

u^Vi + w ;(f 3 f 5 ' + y 3 V 6 ) + ™ 3 (y 3 V 3 + f 5 V 6 



■to, 



-™ 3 Y 4 V 4 + w' 6 (Y 3 Y^ + Y 3 %) + w' A {Y 3 Y 3 + y 6 V 6 ; 

{utoOta - nn') + «* [y;(y - y 6 ) + y 2 (y« - yd] } = 

w' 3 {w 3 (y^y 5 - r 6 r 5 ') + Wb [y 2 '(y 6 - r 5 ) + y s (^ - n')] } 
-«i {tuB^Vs - y 2 y 2 ) + W3 [y 5 (y 8 # - y 2 ) + r 5 '(r 2 - y 3 )} } = 
«i {«*(y 2 >J - Y s Yi) + w 5 [y 5 '(y 3 - f 2 ) + f 6 (f 3 ' - r 2 ')] } 
-«i {^ 6 (f 6 V 6 - r 5 r 5 ') + W4 [f 2 '(f 5 - y 6 ) + y 2 (Yq - r 5 ')] } = 
w' 4 {w 4 (y^y 5 - r 6 r 5 ') + W6 [y 2 '(y 6 - y s ) + y 3 (y^ - r 5 ')] } 
4 {w 6 (Y 3 % - y 2 y 2 ) + w 4 [y 5 (y 3 - r 2 ) + r 5 '(r 2 - r 3 )] } = 
i {w 4 (y 3 % - y 3 y 2 ) + w 6 [y 5 '(y 3 - f 2 ) + r 6 (r 3 ' - r 2 ')] } 



(A.5) 



(A.6) 



(A.7) 



(A.8) 



(A.9) 



(A.10) 



(A.ll) 



The functions Wj(X) are some of the Boltzmann weights of D 2 model and are given by JIG 

q " ^(A) = ~(ff--)(^-4)(e _A + l) 
2 q q 



w 2 (X) = (e A - e- A )(- 



Q 



1 



Is ,e 



9 



l,,e A 



5 



«*(A) = --(<? - -)(- - -^-)(e A + 1), «, 5 (A) = -(q - -)(- - -^-)(-e A + 1) 



q q 



q q 



q q 



w 6 (X) ■ 

Appendix B : A new D 2 .R-matrix 



2 '" (/ (/ e 
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(A.12) 

(A.13) 
(A.14) 



We begin by presenting the new D 2 i?-matrix 



R(X) 



(e 2A - q 2 ) 2 



a^2,3 



A3 



(g 2 -l)(e 2A -g 



£ aa + g(e 2A -l)(e 2A -g 2 ) £ 

"/ft?' 

aor/3 ^2,3 

-[(e A + l) ^ +e A (e A + l) £ }(E a0 ® ^ Q + ^ v ® £ a y ) 

a < 2 a > 3 

= 2,3 = 2,3 
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(g 2 -l)(e 2A -g 2 , r/1 . _ 



[(l-e A ) £ +e A (e A -l) £ ](£ a/3 ® ^ Q + ^ v <g> 

a < 2 a > 3 

= 2,3 /9 = 2,3 

+ E ® ^ + b+(\)E p > a , <g> ^ + &~(A)£ Q/3 <g> £7 aV3 + 6-(A)S^ ® E Pa 

a ^ 2,3 
^ = 2,3 

+ X! c + (A)E aa ® £ aV + c-(X)E aa ® £ aa + d + {\)E aa , ® E o!a + d-{\)E aa , ® £ aa * 

a=2,3 

+ ^ /(A) [5 M - ® E aa + £ aa ® £ ac/ - £ aa ® - E a > a ® 

a=2,3 

(B.l) 

where £7 a/ g are the elementary 4x4 matrices and we set a' = 5 — a. The Boltzmann weights 
are given by 

a n (A) = a 44 (A) = g 2 (e 2A - l) 2 , a 14 (A) = a 41 (A)e~ 2A = (q - l)(g 2 - l)(e 2A + q) (B.2) 



3/2 — 1/2 

6f = ±^-(g 2 - l)(e 2A - l)(e A ± 1), bf = ±^-(g 2 - l)(e 2A - l)(e A ± q 2 ) (B.3) 
6f = ^e A (g 2 - l)(e 2A - l)(e A ± q 2 ), 5f = ^e A (g 2 - l)(e 2A - l)(e A ± 1) (B.4) 
& = - 1) V ± I) 2 , /(A) = ^(e 2A - l)(g 2 - l) 2 (B.5) 

c± = ±^(g 2 - l)(e A =F l)[e A (3 + q 2 ) ± (1 + 3g 2 )] + g(e 2A - l)(e 2A - g 2 ) (B.6) 

This i?-matrix has additional Boltzmann weights, the last term in equation (B.l), as com- 
pare to the standard D\ i?-matrix ]I0| . In addition, several other weights have also a different 



functional dependence on the spectral paramater A. For periodic boundary conditions, such 
differences are not important since we verified, by using the algebraic Bethe ansatz approach 



[ |i~5|1 , that the corresponding Bethe ansatz equations and eigenvalues are the same as those 
found for the standard D\ model |T6|, [OJ. This result is a strong indication that indeed the 



i?-matrix (B.l) can be obtained by twisting the usual D\ i?-matrix. However, the situation for 
open boundary conditions turns out to be a bit different. In fact, we did not find any diagonal 
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solution of the corresponding reflection equation. The basic K-matrices are non-diagonal and 
we managed to find two classes of such solutions. The first family depends only on a discrete 
parameter e — ± and is given by 



n (1,£) (A,£-; 

y 2 (1 - £) (A,£-: 

n (1 ' £) (A,e-; 



^ 3 (1 ' £) (A,£-) = ^(l + ^)(l + e 2A ) 
^ 6 (1 ' £) (A,e-) = i(l-^)(-l + e 2A ) 



(B.7) 
(B.8) 

(B.9) 



while the second family has an extra continuous parameter £_ 



(2,e 



Y, 



(2,e 



Y, 



(2,e 



(2, £ 



(2,e 



(2,e 



} (a,^; 
} (a,£-: 

} (A,e-; 

} (a,^: 

} (A ; e-; 



,2A 



e_)(q + ee 2X )e 



(l-Ce^)(q + ee 2X )e x 
1 



;i-e-e A )(e- + e A )(l + e 2A )(e + g), 



-(l + ^e A )(-£_ + e A )(l + e 2A )(£ + g), 



- e 
2 



2A 



,2A 



l)( £ - g )(£_-e A )(£_e A -l), 
l)( £ -g)(£- + e A )(£_e A + l), 



(B.10) 

(B.H) 
(B.12) 

(B.13) 

(B.14) 

(B.15) 



Finally, we remark that since this new i?-matrix is only unitary, the associated K + (X) ma- 
trices can not be directly obtained by an isomorphism of the type described in (20). However, 
as shown in ref. plj, unitarity is a sufficient condition to allow one to construct commutative 
transfer matrices leading to open spin chains. In this case one has to solve an extra reflection 



equation to obtain the K + (X) matrix EI 
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